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We present a theoretical framework to analyze fluctuations of the electric current through a non- 
collinear single-electron spin- valve transistor in the limit of weak tunnel coupling. The system under 
consideration consists of two tunnel junctions that connect a small, nonmagnetic metallic island to 
two ferromagnetic leads with noncollinear magnetization. We study the current noise spectrum 
as a function of bias voltage, frequency, and the relative angle between the leads' magnetization 
directions and find that both the zero- and the finite-frequency current noise are strongly affected 
by charging energy and spin accumulation in the island. 
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I. INTRODUCTION 

The continuous trend of miniaturization in electronics 
after the realization of the first transistor 1 lead to de- 
vice dimensions that, nowadays, approach the nanome- 
ter scale, at which Coulomb-interaction effects and quan- 
tum mechanics become important. A paradigmatic sys- 
tem for the emergence of such phenomena is the single- 
electron transistor (SET), in which the continuous move- 
ment of charge carriers is replaced by a discrete charg- 
ing and discharging of a small central electrode^ Elec- 
trodes and central island of the SETs maybe composed 
of different materials, involving normal metals, super- 
conductors, and/or ferromagnets. The use of ferromag- 
netic components give rise to spintronic effects such as 
the tunnel-magneto resistance (TMR) that can be used 
in information technology.— — Therefore, SETs based 
on ferromagnetic materials have been extensively stud- 
ied experimentally^ — and theoretically^— Since the 
signal-to-noise ratio is crucial for possible applications it 
is important to investigate current fluctuations. Further- 
more, the current noise can reveal additional informa- 
tion about the transport processes that is not contained 
in the mean current. The signal-to- noise ratio is rep- 
resented by the Fano factor F = S I][ /(2eI), where S 11 
(properly defined below) is the noise of the current / and 
e the elementary charge. Especially in the context of 
mesoscopic devices (such as SETs) current fluctuations 
attracted much interest over the years, for a review see 
Ref. EH- 

It has been shown that interaction effects can de- 
stroy the typical fermionic suppression of the classical 
Poisson noise. Super-Poissonian Fano factors have been 
observed in devices containing electron reservoirs that 
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are coupled to quantum wells, 6 quantum dots 
gle molecules^— carbon nanotubes / 4 ! 77 single-barrier 
semiconductor heterostructure^ and quantum rings* 7 ^ 
In the present work we study current fluctuations of a 
single-electron transistor composed of a central metallic 
island that is weakly tunnel coupled to two ferromag- 
netic leads, whose magnetization directions enclose an 
arbitrary angle 0, see Fig. [TJ This so called noncollinear 
single- electron spin-valve transistor shows typical single- 



electron and spintronic phenomena, such as Coulomb os- 
cillations, Coulomb blockade, TMR, and spin accumula- 
tion. Furthermore, the system exhibits an interaction- 
induced exchange field that exists between the central 
region and the leads i 45 i 46 i 61 This fictitious field is evoked 
by virtual-tunneling processes between the interacting 
central region and the polarized leads and results in a 
precession of the accumulated island spin. The exis- 
tence of such an exchange field in mesoscopic conductors 
was firstly theoretically described^ and experimentally 
confirmed^! - — in the context of quantum dots coupled 
to ferromagnetic reservoirs. 

In literature, there are several publications dealing 
with current fluctuations in special limits of the single- 
electron spin- valve transistor. The complexity of the re- 
quired theoretical description drastically simplifies if one 
considers unpolarized leads (p = 0) only. In this limit 
(two normal leads coupled to a central normal region with 
a continuous level spectrum, NNN) zero- frequency^"— 
as well as finite- frequency^ 7 - — current fluctuations have 
been addressed theoretically. Since in experiments the 
shot noise is often superimposed by other sources of noise, 
e.g. the noise of the used amplifiers or the 1/f noise 
that occurs due to defects in or near the junctions, there 
are just a few experimental works studying the current 
fluctuations of the NNN system. Most of them have to 
restrict their noise measurements to the regime of large 
bias voltages ^2r— However, recently Kafanov and Dels- 
ing measured the noise of the NNN system over a wide 
voltage ranged Most of the published works considering 
current fluctuations of single-electron spin- valve transis- 
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FIG. 1: (Color online) A metallic island is tunnel coupled to 
two adjacent ferromagnetic leads with noncollinear magneti- 
zation directions enclosing an arbitrary angle <j). This system 
is called single-electron spin- valve transistor. 
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tors, i.e., SETs involving ferromagnets, restrict the in- 
vestigations to collinear setups (p ^ 0, <j) G {0,7r})j22ri2S 
The charge- and the spin-current noise has been con- 
sidered for zero as well as for finite frequency. There 
are works that additionally focus on the investigation of 
the effect of spin-flip scattering on the charge- or spin- 
current fluctuations ,28rifiSl It has been predicted that the 
Fano factor strongly depends on the lead polarization, 
the spin-flip scattering strength, and the contact resis- 
tances. The limit of arbitrary angle (j) but absence of 
Coulomb charging effects on the central electrode has 
been studied in Ref. llQll . 

In the present work, we derive a theoretical framework 
relying on a diagrammatic real-time approach that in- 
corporates the general description of noncollinear lead 
magnetization directions and Coulomb charging effects 
on the central electrode. The used theory allows for a 
systematic expansion in the tunnel-coupling strength T. 
Due to the weak coupling between island and leads we 
perform a perturbation expansion of the transport prop- 
erties up to first order in T. The presented theory allows 
for the investigation of the zero-frequency as well as the 
frequency-dependent current noise for the noncollinear 
single-electron spin- valve transistor. 

II. MODEL 

The system under consideration is the single-electron 
spin-valve transistor that is illustrated in Fig. HJ Its 
Hamiltonian takes the form 

H= H r +Hi + H c + Ht '- W 

r=L,R r—L,R 

The left (r = L) and right (r = R) ferromagnetic leads 
are described as reservoirs of noninteracting electrons by 

H r = € rks &\hsv a rksv 7 (2) 
ksv 

with a^ su being the annihilation (creation) operator 
of lead r, momentum fc, and transverse-channel index 
v = 1,2,...,7V C . The majority (minority) spin states 
are quantized along the lead magnetization direction n r 
and denoted by the spin index s = +(— )• The two 
vectors and ur enclose the angle <j>. For simplic- 
ity we choose the density of states p r s , that describes 
the spin-a electrons in lead r, to be energy independent. 
As a consequence, the lead's degree of spin polarization 
p r = (p+ — p r -)/(p r + + p r -) is constant in energy. 

The two contributions Hi and He represent the metal- 
lic island with an energy spectrum s\ that is characterized 
by a typical level spacing As. In our work, we are inter- 
ested in the limit /c^T, eV ^> Ae, where the spectrum 
can be viewed as continuous. The first part 

H i = ^2 e i c Lv c i^i ( 3 ) 

lav 



describes the kinetic energy of the electrons with spin a 
occupying the island level I in the transverse channel v. 
The annihilation (creation) operator of island electrons 

in the state lav is denoted by c^J u . In our model the 
levels of the spectrum are assumed to be independent 
of spin a and channel v. Due to the noncollinear lead 
magnetization directions there is no canonical choice for 
the spin quantization axis of the island. As demonstrated 
in Ref. [6l| it is convenient to choose the island spin- 
quantization axis x\s parallel to the accumulated island 
spin S. In Fig. [2] we introduce two angles that determine 
the orientation of S relative to the lead magnetization 
directions. The angle a is enclosed by the (nz,— n#)-axis 
and the projection of S onto the (n^, n^-plane, and f3 is 
defined as the angle between S and the (n^, n#)-plane. 

The second term describing the island models the 
Coulomb interaction of electrons that occupy the island: 

H c = E C (N - 7V ext ) 2 . (4) 

Here, N is the number of island electrons and the pa- 
rameter Ec = e 2 /(2C$]) represents the charging-energy 
scale of the system with e being the elementary charge. 
The total capacitance = Cl + Cr + Cq is the sum of 
the capacitances of the left and right interfaces and the 
gate. We assume equal junction capacitances (Cx = Cr) 
since they are much less sensitive to the geometry of 
the tunnel contacts than the tunnel couplings. Fur- 
thermore, the bias voltage V is applied symmetrically 
to the leads, i.e., their electrochemical potentials are 
PL = eV/2 and pr = —eV/2. As a consequence, the 
external charge eN ext = CqVg depends only on the 
gate voltage Vq (and not on the bias voltage V). For 
later convenience, we define the energy of the decou- 
pled island system in a microscopic island state \ as 
E x = (x\ (Hi + He) \x) an d additionally we define An 
as the difference of charging energies of N + 1 and N 
electrons, i.e., A N = E C [2(N - A^ ext ) + 1]. 

Electron tunneling between the ferromagnetic lead r 
and the central island is described by the tunneling 




FIG. 2: (Color online) Scheme of the relative orientation of 
the polarization directions of the two leads ul^Hr and the 
accumulated spin on the island S, parametrized by the angles 
a and /3. For clarity reasons, the island spin S is decomposed 
into two parts S = Si n + S pp , where Si n and S pp are the 
contributions in the (nz,, ii#)-plane and perpendicular to it, 
respectively. 
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Hamiltonian 

H T ,r= £ V; a al ksv c lav + -&.c, (5) 

klscrv 

where we already assumed that the tunneling-matrix el- 
ements Vg a are independent of momentum k and trans- 
verse channel index v. Due to the chosen quantization 
axis of the island spin the elements Vg a consist not only 
of the spin-independent tunnel amplitude t r . In fact they 
additionally contain matrix elements of an SU(2) rotation 
that connects the in general different spin quantization 
axes. In terms of the angles a and (3 the tunneling-matrix 
elements of the left lead are given by 
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,(6) 
•(7) 



The elements describing the right lead can directly be ob- 
tained by replacing L — >• R and — >• — 0. The tunneling 
rate for electrons from lead r with spin 5 into the island 
spin state a is given by T r sa /h = 2tt p r s \Vg a \ 2 / h. Finally, 
we define T r a = £ s T r sa , T r = £ a and r = £ r T r . 



III. METHOD 

In Ref. |[l we presented a diagrammatic real-time tech- 
nique to calculate the average current through a non- 
collinear single-electron spin- valve transistor. The ap- 
proach incorporates noncollinearity of the lead magneti- 
zation directions as well as Coulomb charging effects on 
the central electrode that is treated nonperturbatively. In 
order to address the frequency-dependent current fluctu- 
ations of the device, we need to extend this theory ac- 
cordingly, which is described in this section. 

The section is divided into three parts covering the 
derivation of the density matrix that describes the single- 

of the charge- 



electron spin- valve transistor (Sec. IIII A|) 
current (Sec. IIII B|) . and of the frequency-dependent cur- 
rent noise (Sec. IIII C]) . 



A. Reduced Density Matrix 

Since the leads are considered as equilibrium reservoirs 
of non-interacting electrons, we integrate them out and 
derive an effective description, which only contains the 
degrees of freedom of the metallic island characterized 
by the reduced density matrix p re d. For time-translation 
invariant systems its time evolution can be expressed in 
terms of the reduced propagator H(t — to): 



Pred(t) = - to)pred(^o)- 



(8) 



The initial density matrix p^ d of the system is given by 
Pred = limto^-oo p r ed(^o)- In the long-time limit, the sys- 
tem looses any information about its initial state. Hence 



we are free to define the elements of the initial density 
matrix as (#ed)*2 = *xi,xo*X2,xo with Xo being an ar- 
bitrary state of the reduced system. The microscopic 
island states x are determined by \x) = \{^iau}}^ with 
niav = 0,1 representing whether the corresponding island 
level Iv is occupied by a spin-a electron. Furthermore, 
we define the stationary reduced density-matrix elements 
as = (xi| Pred |X2)- In the used notation, the diag- 
onal elements P* = P x correspond to the occupation 
probabilities of state x- They fulfill the normalization 
condition x ^x = I- Eventually, we get the following 
equation for the elements of /) re d : 



pxi 

Xi 



to 



lim U(t 



L °)X2 Xo' 



(9) 



the transform II (u) = 
i0+)t]/h we switch into frequency 



By means of 
J °° dtU(t) exp[— i(uj 
space which is convenient in the following. Then the 
Dyson equation U(uo) = U^ \uj) + (uo)W (uj)U(uj) 
yields the following form of the reduced propagator: 



u(cu) = (n.W(u)) 1 - w 



(10) 



with the free propagator U^°\uj) and the kernel W(cj). 
The matrix elements of the former are given by 



n (0)XiX ,M 
X2X2 v ; 



.x'A, 



ho + i0+ ' 



(11) 



They describe free propagation in time while the elements 
^X2x fl ^ characterize transitions between the matrix el- 
ements P*, 1 and Py*. The kernel W can be calculated 
within a diagrammatic real-time technique that allows 
for a systematic perturbative expansion in the tunnel- 
coupling strength r i 8Q i 1Q2 ~ Ji£ In the present work, it is 
sufficient to expand it up to the first order since we want 
to describe weak coupling between island and leads. In 
this situation, sequential-tunneling processes are domi- 
nant, however, the formalism is generally formulated and 
not restricted to the sequential-tunneling limit. We note 
that in the diagrammatic language, the transformation 
into frequency space incorporates an additional bosonic 
line carrying the energy Huj into the diagrams. 

By applying the final value theorem lim^^o^ + 
+ )II(u;) lim^oo U(t) = p red to Eq. (HO]) we, finally, 
obtain the matrix form of the generalized master equation 
that determines p re d in the stationary limit 







n (0) (cj = o)- 1 -w(cj = o) 



Pred- 



(12) 



Written in component form this equation reads 



dt X2 



X1X2 



Although we have already traced out the lead degrees of 
freedom, due to the continuous island spectrum the sys- 
tem of equations that has to be solved is still highly di- 
mensional. In the rest of this subsection we will follow the 



procedure of Ref. [Hll to strongly reduce this large number 
of degrees of freedom to the relevant ones for electronic 
transport. The first step is to get rid of the off-diagonal 
density-matrix elements on the right-hand side Eq. (fT3|) . 
In our system, the only possibility to change the micro- 
scopic state of the island is tunneling of electrons from 
the leads to the central region or vice versa, i.e., we ne- 
glect intrinsic spin-flip processes in the island due to a 
large spin-flip time scale r s f. Tunneling is described by 
the Hamiltonians Ht,t which conserve charge and trans- 
verse channel v. Hence the elements of the reduced den- 
sity matrix P*, that have to be taken into account fulfill 
the condition ^2 icr ni ay = ^2i a ^ lav . Furthermore, due to 
large number of relaxation channels in metallic islands 
we assume that the electron dwell time Td w is larger than 
the energy-relaxation time r er . This results in the fact 
that coherent superpositions between island states that 
differ in the number of electrons occupying the level I de- 
cay quickly such that J2 a ni av = ^2 a n' l(TV . With that we 
already achieved that the first contribution of the right- 
hand side of Eq. ([T3]) vanishes since spin degeneracy of 



the island spectrum yields E x 



E x' 



for all relevant 



states x an d x' • To get rid of all remaining off-diagonal 
elements of p re d in the second contribution of Eq. ([T3]) 
we neglect quantum corrections to spin quadrupole and 
higher moments, i.e, only spin-dipole moments of the is- 



land are taken into account. This yields n\ 
Eq. (fT3|) simplifies to the kinetic equation 







dt X2 



, and 



(14) 



To describe the electronic structure of the cen- 
tral metallic island it is useful to introduce the spin- 
dependent electrochemical potentials of the island spins 
fjL a . These quantities describe the two spin subsystems 
that are in general out of equilibrium, i.e., /i^ ^ /i^. 
In the considered limit r s f ^> T er ,Td w ? the two differ- 
ent spin species act like two independent reservoirs of 
electrons, which both may be described by the Fermi 
distribution f(E - fi a ) = [exp(^^f ) + l]" 1 . In gen- 
eral, the spin-dependent chemical potentials fi a depend 
on the number N a of spin-a electrons occupying the is- 
land (fi a = fi a (N a )) and have to be determined by the 
expression N a = Yli f \ e i ~ /^Wr)]- But in the case of 
small level splittings Ae <C ksT^eV it is reasonable to 
assume [i G to be independent of N a . 

The simplified master equation Eq. ([T3]) still depends 
on the microscopic island states x that contain all in- 
formation about the individual occupation of each is- 
land level. However, most of this information is irrele- 
vant for the evaluation of the spin dynamics, the average 
current, and the current-current correlation function of 
the single-electron spin- valve transistor. The indepen- 
dent degrees of freedom that have to be taken into ac- 
count are the three components of the accumulated island 
spin S represented by its magnitude S = hpi(ji^ — /i^)/2 
and its spacially orientation that is characterized by 



the angles a and f3 as well as the probabilities to find 
N electrons on the island. The latter are defined as 
P/v = ^2 P x $n,n x with N x being the number of island 
electrons in state x- The equations that enable a calcula- 
tion of each relevant independent degree of freedom are 
provided by the kinetic equations of the charging-state 
projector \N) (N\ = Y2 X \x)(x\$n,n x and the total island 

spin operator S = (ft/2) Ewv w«n/ 
the Pauli spin-matrices vector. One obtains 



cL,<?crcr'QcrV, where a is 



dt 

d ,a 



(15) 



0= d^}= E (X2\$\xi)W£$P x .. (16) 

X1X2X' 

The elements of the kernel W are evaluated with the help 
of the diagrammatic rules explained in detail in Ref. [6jJ. 
By inserting the obtained kernel elements and by addi- 
tionally utilizing the previously discussed model simplifi- 
cations we eventually obtain the kinetic equations of the 
independent degrees of freedom P/v and S. The charging- 
state occupation probabilities are given by 

^P N = n^y] [a+ a (Ajv_i)Pjy_i + a~ a (A N )P N+1 

ra 

-a+ a (A N )P N - a-^An-ijPN] , (17) 

while the time evolution of the accumulated island spin 
can be expressed as a Bloch-like equation 



d(S> 
dt 



fd(S>\ 
V dt A 



rel 



d(S> 
dt 



(18) 



The three contributions describe accumulation, relax- 
ation, and rotation of the island spin. Their explicit 
forms are given in Appendix [A] We note that the 
interaction-induced exchange field that exists between 
the ferromagnetic leads and the metallic island is con- 
tained in the spin-rotation term . 

In both sets of the kinetic equations, Eqs. (fTT|) and 
([T8]h the island rate functions 



af a (E) := ±a° r 



E - (fi r - fi a ] 


) 


exp 


, E—^jlr-jlcr) 

^ k B T 


- 1 



(19) 



PiN c 
2-Kh 



ap- 



with the dimensionless conductance a® 
pear. They describe tunneling of spin-a electrons with 
energy E between lead r and central island. 



B. Charge Current 

The charge current operator I r = —edN r /dt through 
lead r, where N r is the total number of electrons in lead 
r, is obtained from the quantum- mechanical equation of 
motion in Heisenberg picture, 



le r 



[H,N r ] = -j V£ajww+H.c.. (20) 



klsav 
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A comparison of this expression with the tunneling 
Hamiltonian of lead r, see Eq. (|5]), yields that I r can di- 
rectly be obtained from i^r,r by performing the replace- 
ment Vg a — >• —ieVg a /h. The constant complex prefactor 
can easily be incorporated in the used diagrammatic tech- 
nique. Eventually, the stationary charge current through 
lead r is given by 

Ir= £ (W'Xi^* ( 21 ) 

XXlX2 

where the matrix elements of the current transition rates 
W /r are directly obtained by multiplying the correspond- 
ing matrix elements of W with the net transported 
charge from lead r to the island. We note that minus 
signs originating from the complex prefactor and factors 
1/2 that appear due to the definition of the symmetrized 
current have to be taken into account appropriatly. 

Since for the average current charge conservation yields 
I = II = —Ir, we can consider the symmetrized cur- 
rent / = (II — Ir)/2, which will be particularly useful 
for the discussion of the finite-frequency noise, see sub- 
section [TlLCl The symmetrized current can be written in 
matrix notation as 

7 = tr[W / (^ = 0)p red ]. (22) 

By adapting the same procedure that we used in deriving 
the kinetic equation Eq. ([T4]) to the expression of the 
current / one gets 

7 = f E ± K( A n) - «-(Aiv-i)] Pn ■ (23) 

rNa 

The upper (lower) sign has to be taken for the left (right) 
lead. 



C. Current Noise 

Fluctuations of the current are described by the 
current- current correlation function 

S H (t) = </(t)/(0)> + (J(O)J(t)) - 2 (I) 2 . (24) 

Its Fourier transform defines the frequency-dependent 
current noise: 

00 

S h (uj)= J At S n (t)e- iut 

— CO 

OO 

= Jdt[(I(t)I(0)) + (J(0)/(t))] e-^-47r<I> 2 5(M}25) 

— OO 

By definition, S IJ[ (uj) is symmetric in frequency, 
S IT (uj) = S it (—uj). It represents a real physical ob- 
servable. For finite frequencies, the total current flowing 
through the single-electron spin-valve transistor is not 



equal to the symmetrized current / = (II — Ir)/2 that 
occurs in Eq. ([25]) . In general, displacement currents ap- 
pear that have to be taken into account by defining the 
current as I = (ClIl — CrIr) / (Cl + Cr)£^ However, in 
our model we assume symmetric junction capacitances, 
see section [III and hence the general formula corresponds 
to the definition of the symmetrized current. 

In the following, we demonstrate how the calculation of 
the noise S ti (uj) can be realized within the diagrammatic 
technique that was used to derive the kinetic equations 
and the charge current formula. In Ref. [64], Braun et al 
developed a formalism relying on the same diagrammatic 
theory that enables to calculate the frequency-dependent 
current noise of a quantum-dot spin valve, which is a 
single-level quantum dot coupled to two ferromagnetic 
leads. However, due to the large number of levels con- 
tributing to transport this theoretical framework is not 
applicable to our system. Hence we present how we ex- 
tend the theory to describe the current-current fluctua- 
tions of the single-electron spin- valve transistor. 

The frequency uo appearing in Eq. ([25]) is taken into 
account in the diagrammatic language by introducing an 
additional bosonic line carrying the energy hw that con- 
nects the two current vertices that replace two vertices 
originating from the tunneling Hamiltonian^ 

In matrix notation the frequency-dependent current 
noise is given by 

S n (uj) = -47r (I) 2 S(uj) 

A tr [W /J ( 7 ^ red 

+W^( 7 ^)n( 7 ^)Wi( 7 ^)p red ] .(26) 

with three new diagrammatic objects W /j (c<j), W<(cj), 
and (cj). 

In W^(a;) and W>(cj) one current vertex placed on 
the upper or lower Keldysh contour is contacted by the 
bosonic line that enters the diagrams from the left or 
leaves them to the right, respectively. In the diagrams of 
W 7/ (co>) the additional line connects two current vertices. 

The noise formula given in Eq. ([26]) represents the gen- 
eral expression of the current fluctuations of the single- 
electron spin- valve transistor. In the limit of weak island- 
lead coupling, only diagrams containing one tunneling 
line contribute to the kernels, i.e., they are expanded up 
to first order in T. For a systematic and consistent per- 
turbation expansion of the expression for the noise, we 
count the frequency as one order in the tunnel-coupling 
strength, uo ~ T, and expand to first order in this small 
parameter. In this limit the frequency dependence of the 
kernels can be neglected since each correction in uj evokes 
contributions that are at least proportional to T 2 . The 
only frequency dependence we keep is contained in the 
reduced propagator II (a;), see Eq. ([TO]) . This expression 
is treated consistently when only the frequency depen- 
dence of the free propagator is taken into account. As a 
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result, we find 
S n (uj) = -47r (I) 2 S(uj) 

+ I^tr [W 7/ p red 



7 =± 



+W / (n^( 7 a;) -W)- 1 W / p red l(27) 



with the definitions W = W(uj = 0), W 7 = W 7 (u = 
0) = W£(cj = 0), and W /J = W /j (cj = 0). This 
formula for S I][ ((ju) depends on all the elements of the 
reduced density matrix Pg* . An effective description 
that only contains the charge-state occupation probabil- 
ities P/v and the accumulated island spin S as degrees of 
freedom is obtained by executing the same procedure as 
in the derivation of the kinetic equations. This enables 
to remove the coherent superpositions on the right-hand 
side of Eq. ([27]) , i.e., only diagonal matrix elements P x 
enter. In the island-charge-state basis the frequency de- 
pendent noise is then given by 



S 11 (w) = -47r (iy 5(u) + e 1 W 77 P 
1 



-e 

2 



^W 7 (n^( 7 a;) 1 _W)- 1 W / 



■7=± 



P, (28) 



where the vector e is defined by e n = 1 for all N and 
the vector of the island-occupation probabilities P = 
(.., P/v-i, Pn, Pn+i, •••) fulfills the normalization condi- 
tion e T P = 1. The matrix elements of the kernels in 
charge space that have to be plugged in the formula 
Eq. ([28]) are defined as follows: 

Wn,N' =/ K^2 {dN,N>+lttrcr(AN-l) + Sn+I^' Oi~ a (A N ) 
ra 

-S n ,n> [a+ (A N ) - a-^AN^)] } , (29) 

W N,N> =e7r ^2± [(5Ar,Ar/ + ia+ a (AAT-l) 



— 8N+l,N'® ra (AN)] , 

, — 6^7T > 

ra 

+8N+l,N'Oi~ a (AN)] • 



(30) 



(31) 



In Eq. ([30]) the upper/lower sign has to be chosen for the 
kernel contributions representing the current through the 
left /right tunnel junction (r = L/R). 

As mentioned above, in the considered limit the only 
frequency dependence that is contained in the noise is 
represented by the free propagator of the system. In 
charge space it is given by 



n (0) Nl d a 

11 N 2 N'\ UJ ) 



-huo + i0+ 



^N 1 ,N f ^N 2 ,N' 



(32) 



Hence II^ ^(a;) is represented by a diagonal matrix. For 
finite frequencies, the z0 + in Eq. ([32]) drops together with 
the delta function in Eq. ([28]) . 



In the Eqs. ([29|)-(|3l]) one finds that the kernels W, 
W J , and W /J are proportional to the island rate func- 
tions afp , which makes the expression of the current noise 
reliable for frequencies lj < af a . 



IV. RESULTS 

In the present section the results of our calculations 
concerning the current fluctuations of a single-electron 
spin- valve transistor are presented. We separately dis- 
cuss the zero-frequency limit and the finite-frequency 
noise of the system in subsections IIV Al and IIV B] respec- 
tively. For simplicity we assume that the ferromagnetic 
leads are symmetrically polarized pl = Pr = P, which 
can be experimentally realized by using electrodes of the 
same material. However, the tunnel-coupling strengths 
to the leads T r are not restricted to a symmetric setup 
and we define the asymmetry parameter a = Tl/Tr. 



A. Zero- Frequency Current Noise 

We start with the case of collinear polarization of the 
ferromagnets. Due to the symmetric degree of polar- 
ization there is no spin accumulation on the central is- 
land for parallelly aligned lead magnetization directions 
(4> = 0). This results in the fact that both the current 
and the noise are independent of the lead polarization 
and correspond to those of the unpolarized case i 84 i 85 In 
the more general setup Pl 7^ Pr there is a finite island 
spin accumulation which results in a polarization depen- 
dent Fano factor. The situation changes for the antipar- 
allel setup (<fi = tt) as in this case the spin accumulation 
strongly depends on p. In Fig.[3]we show for the antipar- 
allel case (a) the current, (b) the Fano factor, and (c) 
the second derivative of the current for different values 
of the lead polarization and a large asymmetry parame- 
ter (a = 10). We want to emphasize that for obtaining 
reliable results in the Coulomb-blockade regime, where 
sequential-tunneling transport is completely blocked (low 
bias voltages), cotunneling processes have to be taken 
into account. This regime is marked by the grey area. 
For vanishing polarization the sub-Poissonian Fano fac- 
tor exhibits strongly pronounced Coulomb blockade oscil- 
liations which represent the steps in the current-voltage 
characteristics. Hence the maxima of these oscillations 
occur at the excitation energies of the relevant island 
charging states of the central island. The sharp structure 
of the p = Fano factor smears out for finite polariza- 
tions. Additionally, the position of the maxima change. 
Both effects are caused by the island spin accumulation 
since it induces the spin dependence of the charging-state 
excitation energies that results in a more complex exci- 
tation spectrum. To illustrate this behavior we plotted 
the second derivative of the current in Fig. [3fc) and con- 
sider the case of p = 20%. The Coulomb peaks of the 
unpolarized situation are split up into two peaks rep- 
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eV/Ec 

FIG. 3: Antiparallel setup (0 = tt) of the single-electron spin- 
valve transistor: (a) average current, (b) Fano factor, and (b) 
second derivative of the current as a function of the applied 
bias voltage V for different lead polarizations p. For all three 
plots, the remaining parameters were chosen to be CgVg = 
3e/10, E c = 50/c s T, and a = 10. 



resenting the two different spin reservoirs of the island 
(marked by the up and down arrows). We emphasize 
that the first peak (No + 1) does not split due to the 
exponentially suppressed island spin accumulation in the 
Coulomb-blockade regime. For the charging states No+n 
with n £ N + the spin-down (spin-up) excitation ener- 
gies are shifted towards lower (higher) bias voltages and 
vice versa for N — n. However, due to the asymmetric 
coupling the occupation of the latter states is strongly 
suppressed. 

Even for very large asymmetry parameters a ^> 1 and 
highly polarized materials p > 90% the obtained Fano 
factors remain sub-Poissonian. This is in contrast to sys- 
tems consistent of two ferromagnetic leads tunnel cou- 
pled to a central region that exhibits a discrete energy 
spectrum i 64 i 71 i 96 In the case of a single-level quantum 
dot coupled to two ferromagnetic leads (quantum-dot 
spin valve) super-Poissonian statistics arises as a result of 
bunching effects that are caused by spin blockade. In the 
single-electron spin- valve transistor such bunching effects 
do not occur since its transport behavior is significantly 
different. To illustrate this we start by considering the 
quantum-dot spin valve with parallel aligned lead mag- 
netization directions in the regime where the single level 
of the dot is predominant occupied by one electron. Fi- 



nite polarizations evoke that mainly majority spins enter 
the dot. These can easily leave into the drain lead due 
to the parallel lead magnetizations. However, if a mi- 
nority spin occupies the single level then its large dwell 
time (few minority states in the drain electrode) leads to 
a temporary blockade of the current through the system 
since there are no further charging states in the trans- 
port window. Hence the rare event of a minority charge 
carrier tunneling onto the dot bunches the flow of ma- 
jority electrons. This behavior leads to super-Poissonian 
noise. In contrast to this, the continuous level struc- 
ture of the single-electron spin-valve transistor prevents 
bunching caused by minority electrons. This is due to the 
fact that there are two different processes that change the 
charging state of the system after a minority electron tun- 
neled onto the island. On the one hand the same electron 
or another minority charge carrier can leave the central 
electrode into the drain with a small probability (sim- 
ilar to the quantum-dot spin valve). But on the other 
hand, contrary to the single-level system, it is addition- 
ally possible that also majority electrons can tunnel into 
the drain. Even for higher polarizations the latter process 
is not suppressed. It transfers the system into the lower 
island charging state and enables subsequent repetitions 
of island filling and depletion processes. Therefore, the 
transport through the system is not blocked by minor- 
ity electrons as in the case of single- level quantum dots, 
arising in sub-Poissonian Fano factors. 

Having discussed the limit of parallel and antiparallel 
magnetization directions we now turn to the noncollinear 
single-electron spin- valve transistor. If one neglects the 
Coulomb interaction of electrons on the central island 
and additionally assumes that the leads are held at zero 
temperature then the following analytic expression of the 
Fano factor is obtained*^ 

F=i(l+ P 2 sin 2 £). (33) 

In this limit, the Fano factor exclusively depends on the 
polarization of the leads and on the angle <f) between the 
p- vectors. The formula represents a monotonically in- 
creasing F from its minima F(0) = 1/2 to its maximal 
value F(tt) = (1 +p 2 )/2. However, by neglecting the 
Coulomb-repulsion energy crucial single-charging effects 
that govern the transport characteristics of the system 
are not taken into account. Hence Eq. ([33]) is an insuf- 
ficient description of the noise in a single-electron spin- 
valve transistor. 

In our formalism the electron-electron interaction on 
the central island is taken into account nonperturba- 
tively. Its interplay with the finite spin polarization is 
giving rise to the exchange field that acts on the accu- 
mulated spin on the central island. In the following the 
impact of the field on the current fluctuations of the sys- 
tem is discussed in detail. The exchange field is an imma- 
nent part of the used theoretical framework and occurs 
in the kinetic equation of the island spin, see Eqs. (fT8|) 
and (|A3p . The expression of the exchange field contribu- 
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FIG. 4: Current through the noncollinear (0 = 7r/2) 
single-electron spin-valve transistor with (solid) and without 
(dashed) the effect of the exchange field. The dashed verti- 
cal lines represent the threshold voltages of different charging 
states. The remaining parameters were chosen to be p — 9/10, 
CgVg = 3e/10, E c = 50k B T, and a = 1. 



tion that is induced between island and lead r is given 
by Eq. (|A4|) . Due to the lead Fermi functions that ap- 
pear in the integral the exchange field depends on the 
applied bias voltage. Therefore, to investigate how the 
exchange field affects the current noise we have to identify 
a voltage regime where the field noticeably influences the 
transport. To this end, we consider the current through 
the system as a function of the bias voltage for a non- 
collinear angle <p = 7r/2, see Fig. 2J In the presented 
figure the influence of the exchange field is illustrated by 
the difference of the two shown graphs. Here, the dashed 
line represents the current of an artificial situation where 
both exchange-field kinetic-equation contributions B£ xc 
were manually set to zero. While for the calculation of 
the solid line the field is fully taken into account. A 
comparison of the two graphs yields that in the vicinity 
of the threshold voltage, that enables the occupation of 
the charging state No — 1, the field strongly affects the 
current- volt age characteristics. Here, we point out that 
to observe a significant effect of the exchange field highly 
polarized leads with p > 0.7 have to be considered. 

After identifying a voltage regime where the exchange 
field is pronounced we now describe its influence on the 
current fluctuations of the single-electron spin- valve tran- 
sistor. In Fig. [5] the charge current, the zero-frequency 
noise, and the Fano factor are plotted as a function of 
the angle between the lead polarization directions <fi. All 
three plots are calculated for symmetric tunnel-coupling 
strengths (a = 1) and a bias voltage of V = 2Ec/e 
that lies in the voltage window exhibiting pronounced 
exchange fields. The precession of the island spin caused 
by the field leads to an increase of current and noise. 
However, the enhancement of the latter is weaker for all 
angles <fi and therefore the spin precession results in an 
decreased sub-Poissonian Fano factor. This behavior is 
also present in systems with asymmetry parameters a be- 
ing unequal to one. The plotted Fano factors exhibit a 
monotonic behavior between the two collinear situations. 
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FIG. 5: Noncollinear single-electron spin- valve transistor: 
(a) charge current, (b) zero- frequency noise and (c) Fano fac- 
tor as a function of the angle between the lead polarization 
directions cj) with (solid) and without (dashed) the effect of the 
exchange field. The parameters were chosen to be p — 9/10, 
eV = 2E C , CgVg = 3e/10, E c = 50k B T, and a = 1. 



We emphasize that this is caused by the special choice of 
the parameters, i.e., there are parameter sets for which 
the maxima (minima) of F do not occur at the collinear 
angles. In Fig. Efc), we additionally plotted the graph of 
the Fano factor of the noninteracting situation (Ec = 0) 
given by Eq. ([33]) . The comparison shows that a consid- 
eration of the noninteracting limit is insufficient to obtain 
reliable results of the current fluctuations of the single- 
electron spin- valve transistor. Even for collinear setups 
the interacting Fano factor strongly deviates from the 
noninteracting one. 



It is worth to mention, that in the exact parallel and 
antiparallel situations the exchange field contributions of 
both ferromagnetic leads point in the same direction as 
the spin accumulation and hence the spin does not rotate 



in Bf 



B R 



i.e., the respective graphs coincide. 
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B. Finite-Frequency Noise 

To analyze the frequency dependence of the current 
noise we consider F(uj) for different bias voltages. The 
result is shown in Fig. [6] We restrict our investigations 
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FIG. 6: Current fluctuations of the noncollinear (0 = tt/2) 
single-electron spin-valve transistor: (a) Zero-frequency Fano 
factor as a function of bias voltage V. (b) Frequency- 
dependent Fano factor as a function of uj for different V. 
For both plots the parameters were chosen to be p = 2/10, 
CgVg = 3e/10, E c = 50/c^T, and a = 40. 

to a bias voltage range where the island is exclusively 
occupied by No or by iVo + 1 electrons. In the up- 
per plot the zero-frequency Fano factor is shown and 
the different voltages Vi with i = 1,2,3 for which we 
will study the current noise are marked by the dashed 
black lines in Fig. [6]^a). The energy eV\ = 0.84Ec is 
slightly larger than the excitation energy of the charg- 
ing state No + 1, i.e., it represents a Coulomb step in 
the current- volt age characteristics. The other two volt- 
ages are lying in between two charging steps and the 
condition P/v + Pn +i = 1 is fulfilled. The respective 
frequency-dependent Fano factors (fixed Vi) are plotted 
as a function of the frequency uj in Fig. [6fb). For all 
voltages the Fano factors are composed of a dynamical 
frequency-dependent contribution reflecting correlations 
(V2,Vs) and anticorrelations (V±) as well as a constant 
contribution. The explicit value of the latter is given 
by S it (uj ^> Ec/K). A comparison shows that for larger 
voltages the noise spectra are broadened and the constant 
term is shifted below the value of the zero- frequency noise 
to current ratio. 

All three graphs in Fig. [6f b) indicate that the struc- 
ture in the noise is destroyed for large enough frequencies. 
Hence we consider the half width ujf 1/2 °f the frequency 
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FIG. 7: Half width of the frequency-dependent Fano factor as 
a function of V for different polarizations p. The frequency 
is measured in units EcpiN c and the remaining parameters 
were chosen to be <j) — tt/2, CgVg — 3e/10, Ec — 50/csT, 
and a = 40. 



dependent Fano factor as a function of the applied bias 
voltage to determine the scale on which the system looses 
its correlation information, see Fig. [71 In the considered 
voltage regime where only two charging states are al- 
lowed the half width is linear in V . The slope of the 
obtained straight line is proportional to 7rao, with the 
used definition ao = ^2 rcr o^. a • It is maximal for p = 
and deacreases with increasing degree of lead polariza- 
tion. Here, we emphasize that as soon as the voltage is 
large enough to bring additional higher charging states 
into the transport window this simple behavior of oop 1/2 
is not valid anymore. 



V. CONCLUSION 

We have presented a diagrammatic real-time tech- 
nique to investigate the current fluctuations of the single- 
electron spin- valve transistor in the limit of weak tunnel 
coupling. The theory allows for noncollinear lead mag- 
netization directions and simultaneously takes Coulomb 
charging effects on the central electrode into account. 

For the collinear as well as for the noncollinear setup 
sub-Poissonian transport statistics were observed. We 
found that finite lead polarizations smear out the char- 
acteristic Coulomb blockade oscillations of the Fano fac- 
tor as a function of the bias voltage. The origin of this 
behavior was addressed to the splitting of the charging- 
state excitation energies caused by the nonvanishing spin 
accumulation on the metallic island. Furthermore, we 
identified the voltage regime in which the effect of the ex- 
change field is pronounced and demonstrated that B£ xc 
leads to a reduction of the zero- frequency Fano factor. 

Finally, we analyzed the frequency dependence of the 
current noise and found that F(u) reflects correlations 
and anticorrelations. However, this noise structure is lost 
for high enough frequencies. We identified the frequency 
scale for the corresponding crossover. 
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Appendix A: Kinetic equations of island spin 



with the Bohr magneton /ie, the dimensionless mag- 
netic moment of electrons g, and the notation f^~(E) = 
f(E — fi r ) or f~(E) = 1 — f(E — fi r )' Furthermore, 
the definitions of the energy-dependent spin density in 
the island s(uo) = ^[/(uj — fi^) — f{uj — fi^))ris and the 
interaction induced exchange field between central elec- 
trode and lead r 



The three different contributions of the kinetic equa- 
tion for the island spin describing accumulation, relax- 
ation, and rotation processes are found to be 



fd<S)\ 



MS)\ 

v a* J TO ; 



2 L < 1 [r 

Nra 



[^(Aiv-i) - a+ (A N )] P N , (Al) 



= -^2P N duj T r s(uj) 

rel ^ri/ ^ 



[/-(w + Ajv-i) - + A N )] (A2) 



^Jdu sMxB^u), (A3) 



/ \ Pr~Fr-N c a 

=M= o n r 



xjduj' 



N 



uo' — uo 



0/ — LJ ■ 



±N-1 



(A4) 



have been used. A detailed discussion of the different 
spin contributions and the exchange field is presented in 
Ref. EH. 
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